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$_, ■ Abstract 

^\ . A method of a consistent consideration of the phonon contributions to mass and gap operators 

in non-magic nuclei is developed in the so-called g 2 approximation, where g is the low-lying phonon 
creation amplitude. It includes simultaneous accounting for both the usual non-local terms and 
the phonon tadpole ones. The relations which allow the tadpoles to be calculated without any 
new parameters are derived. As an application of the results, the role of the phonon tadpoles in 
the single-particle strength distribution and in the single-particle energies and gap values has been 
considered. Relation to the problem of the surface nature of pairing is discussed. 
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I. INTRODUCTION 



In the last two decades, the progress was achieved in the many-body nuclear theory 
[l|, |2[ in going beyond the standard RPA or the Theory of Finite Fermi Systems (TFFS) by 
means of accounting for coupling of the single-particle degrees of freedom with the low-lying 
collective excitations ("phonons"). For magic nuclei Refs. [3|, |4J, [y, la] could be cited, and for 

111 ] . For the first studies in this field the use 



nuclei with pairing these are Refs. [7|, la, ]9|, ]l0|, 
of the phenomenological Saxon- Woods mean field is typical. In this case, a double set of the 
phenomenological parameters is necessary, the first one for the effective force and the second 
one for the mean field. Note that in Ref. 3| the main terms of the mean field potential well, 
i.e. the ones without the phonon corrections, were found in a self-consistent way, within 
the self-consistent version of the TFFS. However, the phonon characteristics were calculated 
with the use of the Saxon- Woods potential, thereby the double set of the parameters was 
used, too. 

The characteristic feature of the modern developments in this field is the complete refusal 
from the phenomenological mean field when both the field and the phonon characteristics 
are calculated self-consistently, with the use of only one set of parameters for the forces 
to calculate the mean field. The HFB calculations with the Skyrme forces [121 . ]l3( should 
be cited here and the ones, for magic nuclei only, within the relativistic mean field theory 
a, uM- Article [a] is an expressive example. Here, in the framework of the old dynamic 



scheme of [la], the Saxon- Woods mean field was changed by the relativistic mean field. In 
both studies the authors solved the Dyson equation for the mass operator in the so-called 
g 2 approximation, g being the phonon-particle coupling amplitude. Properties of the odd 
nuclei nearby 208 Pb were considered, including the single-particle strength distribution and 
characteristics of the single-particle levels. 

A consistent generalization of the g 2 approximation to the nuclei with pairing has been 
performed in [10]. There Eliashberg's approach [la], which was developed originally for 



superconductivity in the solid state physics, has been used and the realistic calculations 
have been performed for 121 Sn and 123 Sn nuclei. In Refs. 8|, ]l0| the problem of evaluat- 
ing the phonon contribution to the single-particle energies and gap values was consistently 
formulated and the equations obtained were solved. 

The problem of the phonon contribution to the nuclear gap value is, probably, one of 



the most interesting at present because it is in a close relation to the problem of nature 
of the nuclear superfluidity. Note that it is also of great interest for the astrophysics, see, 



e.g., Refs. 17|, |l8j]. As to atomic nuclei themselves, the old problem exists whether the 



nuclear superfluidity has a volume or surface nature. In the latter case, a question remains 
in what extent it appears due to a special form of the initial pairing force or to an additional 
contribution induced by exchange with the collective low-energy surface phonons. The cal- 
culations in [8| have shown that we deal with an intermediate case: the phonon contribution 
to the gap value A for the tin region is about 30% of the experimental gap value, the rest 
of about 70% being due to the "main" pairing force. 

Simultaneously and independently, the problem of nuclear pairing was attacked by the 
Milan group. They combined an ab initio calculation of the gap starting with the Argonne 
V14 NN-force, Refs. [9|, |l9|, with evaluating the phonon contribution to A, Refs. (jj, [9j. Their 
conclusion was that about 50% of the gap value for the tin region appears due to the initial 
NN-force and the other 50% due to the phonon contribution. Note that the alternative 
ab initio calculation of the gap was carried out recently in Refs. 20j (Paris NN-force) and 
2l| (Argonne vi§ NN-force) on the basis of the method developed in [22 ] . Omitting the 
discussion of the two methods of solving the ab initio gap equation, note only that a close 
agreement with the experimental value of A was obtained in Refs. 20j, 21|, leaving a room 
of not more than 20%for the phonon corrections. It agrees qualitatively with the results 
[Si, llQ| , but significantly contradicts those of Refs. |7|, |9( . It is worth mentioning that in the 
cited papers by the Milan group the phonon contribution to the gap of Ca and Ti isotopes 
reaches almost 100%. Thus, a conflicting situation remains in the problem under discussion 
and new investigations are required to clear up the situation. 

However, in all of the above-mentioned papers, and in many other studies dealing with 
evaluation of the low- lying phonon contribution to nuclear characteristics, a term was lost, 
which is of the same g 2 order as the usual non-local pole term of the mass operator taken 
into account. We mean the so-called tadpole term. It was evaluated firstly in the pioneering 
paper by V.A. Khodel 23J and named as the local term. Now we name it, in accordance 
with the particle physics terminology, as the tadpole diagram. This approach is closely 
related to the interpretation of the low-lying surface phonons as members of the Goldstone 
branch which appear due to spontaneous breaking of the translation invariance in nuclei. 
The ghost 1 "-state with the frequency LO\ = is the head of this branch, being the exact 



solution of the self-consistency equation of the TFFS. Therefore, for the simplest version 
of the TFFS effective force, the corresponding eigenfunction is g\ = ^¥-, where U{r) is 
the mean field potential. The TFFS equations for the natural parity excitations will yield 
small excitation energies ul and eigenfunctions gi close to gi, provided the self-consistency 



relation is fulfilled. The main idea of Ref. 23| was to develop such a scheme of evaluating 



the surface phonon corrections to nuclear characteristics that they should vanish for the 
ghost phonon. As it turned out, that is impossible without taking into account the tadpole 
term. 

For magic nuclei the tadpole term was considered in detail in Ref. [2 1 , see also references 
therein. As it turned out, the tadpole contribution to the single-particle energies, splitting of 
the particle- vibration multiplets and other properties of magic nuclei and their odd neighbors 
are, as a rule, important and are often of the opposite sign as compared with the the usual 
non-local terms. Development of an analogous approach for non-magic nuclei is a problem 
of great interest. 

The present paper is the first step in this direction. In Section II we derive a closed 
set of equations for the tadpole operators in nuclei with pairing, which doesn't contain any 
new parameters in addition to those used for calculating the mean-field mass operators. In 
Section III we consider possible approximations to make this set more handy. Closed and 
transparent relations for the particle-hole (ph) and particle-particle (pp) tadpole terms are 
obtained. As an application, modifications of the results of Refs. J8|, llO(] for the single-particle 
strength distribution and single-particle energies and gap values of superfluid nuclei due to 
inclusion of the tadpoles are carried out. The relevance of the tadpole terms to the problem 
of the nature of nuclear pairing is discussed. 

II. PHONON CORRECTIONS TO THE MASS AND GAP OPERATORS 

In this paper, we deal with the nuclei with weak phonon-particle coupling when the 
phonon admixture to the one-particle degrees of freedom could be accounted for within a 
perturbation theory scheme. More specifically, the quantity 

a=- ^-^ (1) 

(2j + IK 



plays the role of the perturbation parameter, and it should be small. Here gi is the average 
value of the matrix element of the L-phonon creation amplitude at the Fermi surface, ujl 
is its excitation energy, and j is a typical value of the single-particle angular momentum 
(j oc A 1 ^). In other words, the coupling strength should be not too high and the excitation 
energy not too low. We call this as the g 2 approximation for mass operators. 

Such a situation takes place in magic and semi-magic nuclei. There is no pairing at all 
in the first case and partially, in the magic subsystem, in the second one. For the magic 
nuclei, the problem under consideration was consistently solved in 3j, see also references 
therein. Our aim is to develop a similar approach for semi-magic nuclei, with pairing in the 
non- magic subsystem taken into account. 

In the general case of nuclei with pairing it is necessary to consider four one-particle 
generalized Green functions, i.e. two Green functions, G and G h , and two Gor'kov functions 
F^ and F^ 2 \ In addition to two normal mass operators, £ and E h , their two anomalous 
counterparts appear. In the textbook [24J thev are denoted as £02 and S 2 o, in Refs. [8|, ll(| 
as E*- 1 ) and Y,^. Here we use the notation [1], where they are denoted as A^ and A^ 2 ). 
Therefore we will often name them as the gap operators. Sometimes we use the term "mass 
operators" for all the four quantities under discussion. 

The main part of the mass operators is determined by the mean field contribution. In the 
g 2 approximation, for the single-particle energy E\ and the average gap value A A we have: 

e x = E AA = ef + <^£ AA (e A ), A A = A AA = A$5 + <^A AA (e A ), (2) 

with the obvious notation. In general, the main, mean field, parts of E and A, in the 
approach discussed, are supposed to be calculated within a self-consistent method, say, the 
HFB one or the self-consistent TFFS. However, in practice the phenomeno logical Woods- 
Saxon potential and pairing force are often used. It should be emphasized that in the 
latter case it is necessary to use the so-called "refinement" procedure [5|, [l0(] in order to 
avoid a double counting of phonon contribution into the terms e\ and the gap A A and 
to extract the "refined" values from the phenomenological ones. In principle, analogous 
precautions should be made in the case of the self-consistent calculation with the use of 
phenomenological forces, too. Indeed, the force parameters should be chosen in such a way 
that the total expressions ([2]) (not the "zero" ones) would reproduce the experimental values. 
Note that this very idea was utilized in 




FIG. 1: Phonon g 2 corrections to the mass operator E(e) in non- magic nuclei (general case) 




FIG. 2: Phonon g 2 corrections to the gap operator A*- '(e) in non-magic nuclei (general case). 

With the use of the g 2 approximation one can write down the phonon corrections, e.g., 
to the operators S and A^\ see Fig.l and Fig. 2, as follows: 



6®X(s) = M(e) + K ph , (3) 

S {2) A (1) (e) = M (1) (e) + K pp . (4) 

The set of diagrams and the relation for S^A^ are absolutely similar to Fig.2 and Eq. (J3J). 
In Fig. 1 and Fig. 2 the empty circles denote the phonon creation amplitudes (vertexes). 
To distinguish between the g- and d- vertexes it is necessary just to look at the direction of 
ingoing and outgoing arrows. In the case of one ingoing and one outgoing arrow we deal 
with the g- (or g h -) vertex. If there are two ingoing arrows, we deal with the S 1 ^ -vertex, two 
outgoing arrows mean the S 2 '-vertex. The terms with two phonon creation amplitudes are 
the energy dependent non-local operators, the sums of them being denoted as M and M^ . 
The last terms are the corresponding ph- and pp-phonon tadpoles. The essential property 
of these tadpoles is that they do not depend on the energy e. 



The explicit expressions for the phonon tadpoles K ph and K pp are given by 

K ph = I ^D(co)Sg(u) , (5) 

K pp = /^HM (1) H, (6) 

J 2"Kl 

where D is the phonon Green function and 5g and 5S- 1 ' are the changes of the ph- and pp- 
phonon creation amplitudes in the external field of another phonon with the same quantum 
numbers as the one whose contribution we analyze. 

A. Magic nuclei 

n 

To begin with, let us first outline briefly the method for magic nuclei, following to [3j. In 
this case there is no pairing and only the first and the last terms in Fig. 1 for the corrections 
to the mass operator E survive. The first term is the usual pole diagram, where the Green 
function G, of course, does not contain pairing effects. The last term means the sum of all 
the irreducible diagrams. In the problem under consideration, as it was mentioned in the 

r i 

Introduction, this sum was evaluated firstly in Ref. [23]. We name it, in accordance with 
the particle physics terminology, as the tadpole diagram. 

The second order in gi correction to the mass operator E reads 

5 { ^(e) = J^S LL (e,u;)D L (co), (7) 

where Dl stands for the L-phonon D-function and Sll is the phonon-particle scattering 
amplitude. As usual, the symbolic multiplication means the integration over intermediate 
coordinates and summation over the spin variables. In accordance with Fig.l, Sll is the 
sum 

Sll = 9hGg L + h9L , (8) 

where 5l9l is the tadpole term. Obviously, it is of the second order correction to E, as far 
as the first order correction to E is the particle-phonon interaction amplitude, gi = 5l^- 
According to the recipe of 3], it can be found by the direct variation of the equation for the 
vertex gL, 

g L = TAg L , (9) 



where T is the Landau-Migdal effective NN-interaction amplitude lj and A = GG is the 
particle- hole propagator. After the variation of Eq. (Q over the phonon creation amplitude 
one obtains: 

6l9l = (8 L F)Ag L + f(5 L A)g L + TA 5 L g L . (10) 

This is an integral equation for the quantity 8lQl with the inhomogeneous term 

m LL = {b L T)Ag L + T{8 L A)g L . (11) 

The procedure of finding the second term of the inhomogeneous term is quite obvious. 
The direct variation of the particle-hole propagator yields: 

(S L A) = 2G(5 L G) = 2GGg L G. (12) 

Up to now, we have used a symbolic notation. To obtain the explicit relations, let us, for 
simplicity, suppose that the mass operator £ is momentum independent. In this case, we 
have for the first order correction 

6%lv{r) = g LM (r) = g L (r)Y LM (n). (13) 

Note that within the Bohr-Mottelson (BM) liquid drop model one has 

,BM (r) = aL <EM_ (u) 

where U(r) stands for the nuclear mean field potential and a^ is the constant which de- 
termines the amplitude of the surface vibration. As it was demonstrated in [3J], the direct 
solution of the RPA-like equation for a low-lying excitation in an even-even nucleus is 
very close to the BM model prescription ffT^|) . Therefore, for qualitative estimations, one 
can use this simplified form of the particle-phonon vertex. 

To obtain the second term in ffTT]) we could fold the quantity ffT2l) with g L . Let us 
introduce the notation T = S^Agi. In the explicit form, with the help of (TT2|) . we obtain 

f de 
T L MiLM 2 { r ,u) = / -—dr 1 dr 2 G{r,r 1 ] e)g LMl {ri) 
J 2txi 

x [G(ri,r 2 ;e-w) + G(r 1 ,r 2 ;e + uj)]g L M 2 (r 2 )G(r 2 ,r;e). (15) 

The problem of finding the first term of f fTTl) looks less obvious. To find it, an ansatz was 

n 

used in [3[ based on the density dependence of the Landau-Migdal amplitude JF, 

b L T = 6 -fs L p, (16) 

dp 
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FIG. 3: Equation for the tadpole in magic nuclei 

where b^p is the transition density associated with the L-phonon excitation. It obeys the 
relation 

S L p = Ag L . (17) 

In the approximation similar to (114p . we have 

(8 L p)™(r) = a L ^-. (18) 

Thus, the equation for the ph-tadpole in magic nuclei reads 

K ph = 5 L TAgD + F(5 L A)g L D + TAK vh . (19) 

In the graphic form this equation is shown in Fig. 3. 

B. Non-magic nuclei 

As it was discussed above, in the general case of nuclei with pairing, in addition to the 
one-particle Green function G, two Gor'kov functions F^- 1 ' 2 ' enter the TFFS relations, and 
two gap functions A^ 1 ' 2 ) appear in addition to the mass operator E. They are related to 
each others, 

A (l,2) = ^ F (l,2) ) (2Q) 

in terms of the interaction amplitude T^ irreducible in the particle-particle channel. It 
should be noted that hereafter we mean that the Green function G takes the superfluidity 
effects into account. 

In the systems with pairing, Eq. ([9]) for the phonon-particle vertex is generalized and two 



new amplitudes appear 



Q: 



dS?00 = M M W = d^\r)Y LM {n). (21) 



Let us omit for a while the subscripts and upper indices. As far as the amplitude T*- in the 



TFFS is considered to be density dependent [26j, [271, |28| , two terms in (121]) appear 



d = J*5F + 6J* F. 



(22) 



It is worth pointing out that usually the first term in (122]) is only taken into account. V.A. 



Khodel 



29] was, evidently, the first who turned attention to the second one. It is natural 



to use for it the ansatz analogous to (TT6"j) : 



Sl^ 






hp, 



(23) 



but now, due to pairing effects, the transition density obeys the equation which is more 
complicated than Eq. (IT7|) . It will be written down below. 

To introduce the standard TFFS notation, let us omit for a while the second term in 
Eq. (|22p . As far as we deal with the low-lying excitations of natural parity, contributions of 
spin-dependent forces could be neglected in the equations for g and g h [l|. As the result, the 
relation g h = g is valid. In this case, the quantities gi and S- 1 ' 2 ' obey the set of equations 
PJ, which could be written in the form similar to (jSJ), 



g = FAg, 



(24) 



but now all the ingredients of (I24p are matrices: 



< 9^ 



T 



\ di2) ; 



'^"00^ 



0^ 
f { 



V 



/ 



A 



C M^ M® 



\ 



\ 



J 



(25) 



(26) 



(27) 



O AA 1 ) m® 
O A?(i) A?(2) 

Here C, M.^, and so on, denote integrals over e of different double products of the Green 
function G(e) and Gor'kov functions, F^(e) and F( 2 \s). They could be found in [lj and 
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we write down now explicitly only two of them, 

C(u) = j ^- [G{e)G{e + u) - F^\e)F {2 \e + w)] , (28) 



and 



Let us come back to the second term of ( 123]) . With the help of the above short notation 
the transition density could be written in a compact form similar to ( TT7j) : 



O(o>) = - J ^- [G{e)F{e + w ) + F(s)G(-e - w )] . (29) 



6p=^A li g i . (30) 



Using this relation, we may express the term under consideration in terms of the "generalized 
vertex function" g: 

M* = ^pX>«ft. (31) 

i 

By substituting this relation to (|23l) we find that the general structure of Eq. (j2"4"|) remains 
valid, but now the "interaction matrix" T becomes more complicated, in particular, non- 
diagonal. To be more exact, the first line of (1261) remains unchanged, but new non-diagonal 
terms appear in two other lines. To simplify their explicit form, let us use the approximation 
for the effective pairing interaction amplitude JF^ which is usually utilized in the TFFS 
(e.g., see |28|). Namely, it is considered as an energy independent delta- force with a density 



dependent strength jF^(p(r)). In this case, the second term of ff22|) is reduced to 

6J*F=^6p(r) X (r), (32) 

where 

X (r)= fp,F(e,r,r) (33) 



2tu 

is the anomalous density. Combining the above relations, one can readily find two new 
non-diagonal terms of the matrix !F, 

^2i = ^3i = ^x(r). (34) 

Thus, we obtain the matrix equation (|24|) for g in the general case where both terms 
of Eq. (122]) are taken into account. After variation of this equation over the field of the 
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surface L-phonon under consideration we find the matrix equation for the tadpole term in 
a superfluid nucleus: 

Sl9l = (hf)Ag L + F(5 L A)g L + TA S L g L . (35) 

In principle, this set of equations solves the problem of finding the tadpole terms under 
discussion. In accordance with Eqs. dSD,®, they should be obtained by folding the solutions 
of Eq. (|35[) with the phonon D-function. However, the explicit form of this equation is quite 
cumbersome. The second term on the r.h.s. of Eq. (|35l) is the most complicated. To find 5£ 
and variations of other components of the A matrix, one can use the well known expressions 
for SG, SG h , 5F^' 2 ^ [1]. In the result, one obtains a lot of integrals of triple combinations of 
the Green and Gor'kov functions of the type of Eq. (TT5|) . To obtain more handy relations, 
some approximations should be made. 

III. SMALL d APPROXIMATION 

As it was discussed in the Introduction, the collectivity of the low-lying ph-phonons, i.e. 
the surface vibrations, exceeds significantly that of the pp-phonons, i.e. the pairing vibra- 
tions. Therefore, we concentrate here on the contributions to the mass and gap operators 
of the phonons of the first type. In this case, the g^ component of the generalized vertex 
(}l dominates in Eqs. (1241 and fl35l) and the inequality g 3> S- 1,2 ^ is valid. For this reason, 
we can omit the terms with the pairing phonon creation amplitudes S 1,2 ^ in these equations 
and, correspondingly, in Figs. 1,2. In this approximation, the diagrams depicted in Figs. 4,5 
should only be taken into account . In addition, we assume that we deal with the "developed 
pairing" case when pairing properties of neighboring even-even nuclei should be considered 
identical. In this case, we have A^ = A^ 2 ) = A and S l) {u) = ±S 2 \-uj) = d(u) 1]. The 
sign "+" takes place for the states of natural parity we consider. The appropriate explicit 
expressions for M(e)=M h (— e) and M^ (e)=M^ (e) are given in [lofl . 

In the approximation under consideration, instead of the set (1241) . we obtain the closed 
equation for g, 

g = TLg , (36) 

and the closed expression for the d- vertex in terms of g: 

d(uj) = (F 2 iC(uj) + J*0(uj) g. (37) 
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FIG. 4: Phonon g 2 corrections to the mass operator S(e) in non-magic nuclei in the small d ap- 
proximation (the Green function G contains pairing effects). 
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FIG. 5: Phonon g 2 corrections to the gap operator A^ l '(e) in non-magic nuclei in the small d 
approximation . 

Sometimes the initial equation (122]) for the d- vertex is more convenient. In the small d 
approximation, it reads: 

d(u) = J*0(u)g + 8J*F . (38) 

The set of Eqs. (135]) for the tadpole terms is also simplified. It could be obtained either 
by omitting terms containing d^ , d^ in (l35j) or by the direct variation of Eqs. (136]) and 
( 1371) . In the result, one receives: 

8g = bTLg + FSCg + TLbg, (39) 

5d = 5F 21 £g + T 2X bLg + 6J*Og + J*dOg + (JF 21 £ + J*0) 5g . (40) 

Thus, we obtained the integral equation (1391) for bg with the inhomogeneous term, which is 
similar to the rrin term (TPT1) for magic nuclei, and expression (HO"]) for Sd. The latter contains 
Sg and four terms, which are analogous to the inhomogeneous terms of the equation for the 
5g. Let us consider them in detail. 

An alternative relation for Sd could be found by variation of Eq. ( |38l) . It is as follows: 



5d = bT^Og + J*&Og + 6J*8F + (8 2 ^)F + T^08g . 



(41) 



To obtain the quantities 8C and SO in the equations discussed above in the explicit form, 
one should variate the propagators C and O: 



S£ = 8(GG-FWfW) 



(42) 
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SO = 5(GF^ +F^G h ), (43) 

and use the small d approximation, omitting the terms with d^ , d^ in the general expres- 
sions for the variation of the Green functions [lj . One obtains: 

SG = GgG - F {1) gF {2) , 5G h = G h gG h - F {2) gF {1) , (44) 

5F^ = GgF® + F^g h G h , 8F&* = F&*gG + G h g h F® . (45) 

For a time, we come back to the notation F^\ F^ 2 \ G h to avoid an explicit specification of 



the energy variables in the integrals similar to that in Eq. IjToj) . which appear after folding 
expressions (1421) and (1431) with g. They could be obtained by combining Eqs. (l42]) - (14~5]) and, 
in the symbolic form, are as follows: 

5Cg = g(GGG - F^F^G + GGG - GF^F® - 

GF (D F (2) _ F W GF m _ F (D F (2) G _ F WG h F^)g , (46) 

and 

- bOg = g(GGF {l) - F (1) F (2) F (1) + GGF {1) + GF {1) G h + 

GF^G h + F^G h G h + F^G h G h - F^F^F^)g . (47) 

Thus, even for the simplified case under consideration, we obtained eight terms for bLg 
instead of the one in Eq. ( TTUl) . In addition to them, eight new terms for SOg appear in the 
expression for Sd. The explicit form of each integral entering Eqs. (146p . (147l) is similar to that 
of Eq. (JI5J) . 

A. Final relationships for the tadpoles 



For magic nuclei, equation ( 1191) for the tadpole term was solved in the coordinate rep- 
resentation [3j. Even in this case the procedure turned out to be quite cumbersome. In 
principle, this method could be generalized to the systems with pairing, using the coordi- 



nate representation for the Green and Gor'kov functions [30|. However, as it is clear from the 
above formulas, in this case it will be much more complicated. For this reason, we prefer to 
use the representation of the single-particle wave functions, the so-called A-representation. 
To make the final equations for the tadpoles more transparent, we also use the diagonal in 

14 




+ ■*- 




FIG. 6: Equation for the tadpole K ph in non-magic nuclei in the small d approximation. 

A approximation. The matter is that the set {A} is chosen in such a way that the mean field 
mass operator £(°) and the corresponding Green function G^ are diagonal in A. We use 
the approximation supposing that the mean field gap function A^°\ the Green function G 
with pairing and Gor'kov functions F^ 1 ' 2 ' are also diagonal in A [lj. I* 1 this approximation, 
C, O and other two-particle propagators contain two A-subscripts, £aiA 2 an d so on, whereas 
in the general case we have £aiA 2 ~~ ¥ £aiA 2 a 3 a 4! etc. The corresponding generalization of the 
equations written down below is quite obvious. 

The equations for the tadpoles are obtained by substitution of Eqs. (1391) and ( HOI) (or 
(HI]) ) into Eqs. ([SJ) and (J6j). The final equation for the K ph tadpole, in the obvious short 
notation, has the form: 

K& = Y] tt-tt- ^1234 (u))£ u (e,u)g M D L (u) + 
*-^ I 2m 2m 

3,4 J 

/ft W fit 1 
2m 2m 

^rwuf^CufrutiK*. (48) 

3,4 J 

For the K pp tadpole, let us first use Eq. (HT1) for 8d. We find: 

K ll = 2 E / ^ L ^ L ^n3A^)0 34 (e,oj)g 3A D L (u;) + 

3,4 J n 

E^u f ^(S0 9 U(e,u)D L (co) + 

3,4 J 

*— ' / Z7TZ Z7U 
3 ^ 

E^^/^^^^X- (49) 

3,4 ^ ^ 

The factor 2 in the first term in Eq. (HOI) appears due to the fact that, in the small d 
approximation, the terms bT^Og and ST^SF 1 in Eq. ( HT1) are equal to each other. 
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FIG. 7: Expression for the tadpole K pp in non-magic nuclei in the small d approximation. 

In the graphic form, Eqs. (HH1) and (|49|) are illustrated in Fig. 6 and Fig. 7 respectively. 
For the sake of simplicity, contrary to Fig. 3, here we do not draw all the internal Green 
functions and omit arrows for those which are drawn. The arrows are only conserved in the 
cases where it is necessary for understanding. In particular, in the last diagrams of both 
figures the arrows show that we deal with the tadpole K ph . Remember that the second 
term on the r.h.s. of Eq. (HHj) in the detailed presentation includes 8 particular diagrams, in 
accordance with Eq. (HEl) . In the case of K pp tadpole, the number of diagrams is even larger. 
Therefore the detailed diagram representation of Eqs. ( 1481) and ()49l) is rather complicated. 

Eq. ( 1491) is convenient for the graphical representation, but it possesses one drawback: it 
contains the term with d^'J 7 ^ with the "hidden" tadpole K ph . To separate the latter one 
explicitly it is necessary to use Eq. PU|) instead of (14T|) . In the result, we find: 

de duj 



•"■12 



E 

3,4 



E 

3,4 



2iri 2ivi 

de du 
2ni 2iri 



<5^1234(w)C , 34(£,^)5 , 34-Dl(w) + 

^(^12)1234 / 2?r . 2?r . (SCg) u (e, uj)D l (u) + 
E-^34 i ^-^-{50gU{E,u)D L {u) + 



3,1 



3.1 



2tti 2ni 



/de f de 

— -£ 3 4(e, UJl) +-7*1234 / -^—.Om{£,wl) 



K: 



ph 

34 



(50) 



One remark should be made concerning the integrals over u in the above equations for 
the tadpoles. The poles of the D-function should be taken into account only because they 
lead to the terms which strongly depend on the low-laying phonon frequency u>l and other 
phonon characteristics. They could change significantly from one nucleus to another. On the 
other hand, the terms appearing due to poles of £, O and other two-particle propagators 
do not practically depend on oul- They are smooth functions of all the variables under 
consideration and should be included into the corresponding mean-field quantities. 
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The solution of the integral equations (HB1) and (|49|) (or (1501) ) yields the tadpole values 
which, according to Eqs. (2) and (3), should be added to the usual non-local terms. Note 
that in the approach discussed there is no need for any new parameters. Below we consider 
some applications of the results obtained. 

B. Applications to description of the single-particle characteristics of non-magic 
nuclei 

In Refs. [8|, [l0| the approach to describe the the single-particle strength distribution for 
non-magic odd nuclei and to take into account the phonon contributions to the single-particle 
energies and gap values has been developed on the basis of generalization of the Eliashberg 



;heory 16] to nuclei, with the first application of the Eliashberg theory to nuclei made in 
3 if . The general set of equations of Refs. [g, [l0| for the energy and gap operators, with 
account for the dynamic spread of a single-particle level (the "dynamic" case), has been 
derived in the diagonal approximation for the mass and gap operators. (Arguments in favor 
of such an approximation could be found in [8|, llO|). The equations are as follows : 



£\r. 



A 



1 + qx v (E Xv ) 

A{ 0) + M^ (E Xv ) 

1 + QXrj(E\ v ) 



Xr, ~ 

Ex, = y/%, + Al, , (51) 

where 

= _Mq^ (52) 

Here M even and M odd are even and odd in energy components of the non-local mass 
operator M (M = M even + M odd ), which enters the r.h.s. of Eq. ©, and M^> is the same 
as in Eq. (J3J). The subscript r\ numerates solutions of the set of Eqs. (l5Tj) . (l52]) . This yields 
the distribution of the single-particle strength in non- magic nuclei. 

In order to obtain the single-particle energies and gap values (the "static" case) from 
Eqs. (i5T|) . (|52|) . it is necessary, for each A, to separate the dominant solution r\ from the set 
{Xtj}. For this aim, the spectroscopic factors should be analyzed. They are given [l0(] with 



,± _ (1 + q\ v )(E Xv ± Sx-n) 
Qx(E Xv ) 



b Xv = A ,„ x 7 I 53 ) 
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where 

& x (e) = (e- 4 0) - M x (e))(e + 4 0) + MJJ(e)) - (Af + M«(e)) 2 . (54) 

The 77-component with the maximal spectroscopic factor should be associated with the 
experimental single-particle level, the details see in [lOJ ] . Let us denote the observed single- 
particle energies and gap values as e\ and A A and the corresponding mean field values as 
e x and A^ . They are related to each other by Eqs. fl51|) . fl52|) with r\ equal to the dominant 
value. Let us rewrite them explicitly, omitting the subscript rj: 

4 0) + Mf en (£ A ) 



£\ 



A) 



l + qx(E x ) 

Af + M} 1 ^ 

l + qx(E x ) 



Ex = \e\ + A* , (55) 



.(o) 



* = -^. (56) 

&X 

are reckoned from the corresponding chemical potentials \x and /tr ' . 



where 

The energies e\ anc 

Note that in Refs. g, [l0|] the phenomenological Saxon- Woods potential was utilized as the 
mean field one and the phenomenological pairing forces were used as well. As far as both 
of them are adjusted to the observed values of E\ and A A , a special "refinement" procedure 
mentioned above is necessary to find e x and A^ values. It is described in detail in the 
cited articles. 

Now it is necessary to modify these results in order to include the tadpoles in accordance 
with Eqs. (151), (fl"]). In fact, there was no specialization of the mass operators in [8|, [lOj to 
derive the relations floTl) and fl55|) . For this reason, in order to include the tadpoles into 
consideration we should just change the mass and gap operators of Refs. [g, ll0[] to the ones 
from Eqs. ([3]) and (jlj). Remember that the tadpole terms K ph and K pp do not depend on 
the energy. Supposing, just as the non-local operators M, M^\ that they are diagonal in A, 
we obtain 

4 0) + M^(E Xv ) t Kf 



Xv 1 + qxr,(E\ v ) 1 + qxr,(E Xr) ) 

Xv 1 + qx v (E\ v ) 1 + qx v (E\ v ) 



Ex, = ^e% + A% , (57) 
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with 

qx v = ^ , (58) 

instead of Eqs. fl5Tl) . fl52l) . In the same way, instead of Eqs. fl55|) . fl56l) for the single-particle 
and gap values, we find: 

£ f + Mr n (E x ) t Kf 

£ X = 7—, 7TT~\ r 



l + q x (E x ) l + q x (E x ) 

X l + qx(E x ) l + q x (E x ) 



Ex = y/% + A\ , (59) 

with 

tt = -^. (60) 

&x 

We see that both the single-particle energy and gap values are changed due to inclusion 
of the tadpoles, both in the dynamic and static cases. In the latter case, the solution of the 
set of Eqs. (1591) should answer the question about the total phonon contribution, including 
the tadpole terms, to the pairing gap, as compared with the mean field, or "refined", value 
A^ . Up to now, all calculations of the phonon corrections to the gap have ignored the 
tadpole contributions. 

Let us briefly discuss the situation in nuclei without pairing. In this case, the equations 
of Sect. II A for the phonon corrections to the single-particle energies were solved in the 



coordinate representation in [3( (see references therein, in particular, [32J). It turned out 
that the tadpole contribution to e x was, as a rule, significant and comparable with that of 
the non-local term of the mass operator. For a qualitative analysis, we limit ourselves with 
the diagonal in A approximation, using Eqs. (|5"T|) and (1591) without any pairing contribution. 
For the spread of a single-particle state we have: 

EXv = ef + M x (s Xv )+Kl h , (61) 

and for the single-particle energies: 

e x = ef ) + M x (e x )+Kl h . (62) 

As far as K ph does not depend on the energy, the shift of the solutions, for a fixed A, will be 
the same for all the values of 77. For the same reason, that is independence of K ph on energy, 
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the spectroscopic factors which are determined by the residues of the Green function and, 
therefore, by the energy derivative of the mass operator in Eq. ([3]), are not changed. 

This conclusion about the role of the ph-tadpole in magic nuclei agrees with the results 



» 



of calculations for single-particle level properties in odd neighbors of 208 Pb in Ref. |6[] cited 
above, where the tadpole contribution was not taken into account. Indeed, as it can be 
seen from Table III of [6j, the authors obtained a good agreement with the experiments for 
the spectroscopic factors, where there is no tadpole contribution, but the agreement for the 
single-particle energies is considerably worse due to the fact that in this case the tadpole 
contribution does exist. 

Things are different in nuclei with pairing. In this case, in the absence of the tadpole, 
the single-particle spectroscopic factors are given with Eqs. (j55|) . (j53|) . If the tadpole terms 
K ph and K pp are included, as it can be easily checked, in addition to a change of E\ v and 
exrj, the expression (1531) itself is modified. Thus, for non- magic nuclei both the energies and 
spectroscopic factors should be changed due to inclusion of the tadpoles. 

IV. CONCLUSION 

In this work, a consistent approach is developed to include the phonon coupling in the g 2 
approximation for mass and gap operators in non-magic nuclei with explicit consideration of 
the tadpole terms, in addition to the usual non-local terms. The general set of equations for 
the phonon corrections under discussion is obtained, which doesn't include any new param- 
eters besides those used in the self-consistent calculation of the "zero" mean field, i.e. the 
ones without phonon contributions. This set is simplified for the case of corrections induced 
by low- lying surface phonons in the "small d approximation" (dP~ ,2 > « g). This approxi- 
mation means that the admixture of pp-phonons with ph-phonons under consideration, i.e. 
the contribution of the d( 1,2 )-vertexes in comparison with the g-vertex, could be neglected . 
The closed integral equation for the ph-tadpole K ph as well as the integral relation for the 
pp-tadpole K pp in terms of K ph are obtained. Even for such a simplified case the relations 
obtained are much more complicated than those for magic nuclei. 

As an application of the relations obtained, the role of the phonon tadpoles in single- 
particle strength distribution, in the single-particle energies and gap values is analyzed. 



Relations of Refs. 



101 ]. where only usual non-local mass operators (in the g 2 approxima- 
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tion) have been taken into account, are modified with the explicit inclusion of the tadpole 
terms. The set of equations obtained is analyzed. Even before numerical calculations, the 
analysis of the structure of these equations and their comparison with those for magic nu- 
clei, lead us to a conclusion that the tadpole terms should change significantly the nuclear 
characteristics under consideration. Indeed, on the one hand, this comparison shows that 
the ph-tadpole K ph in non-magic nuclei should be close to that in magic ones. On the other 
hand, the expression for the pp-tadpole K pp in terms of K ph shows that the first one has no 
smallness in comparison with the second one. Therefore, we could rely on the experience of 
calculations in Ref. [3j for magic nuclei where the contribution of the tadpole term, e.g., to 
the single-particle energies is significant. Note that, contrary to magic nuclei, in non- magic 
ones the tadpoles should also change the spectroscopic factors. A preliminary analysis of 
the modified gap equation shows that here the tadpole could be significant, too. This is 
important for the problem of pairing nature in finite nuclei. 
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